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We consider corrections to the Einstein-Hilbert action which contain both higher order and non- 
local terms. We derive an effective Newtonian gravitational constant applicable at the weak field 
limit and use the primordial nucleosynthesis (BBN) bound and the local gravity constraints on G e / / 
in order to test the viability of several cases of our general Lagrangian. We will also provide a BBN 
constrain on the CR gravitational correction. 
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I. INTRODUCTION 

One would naturally expect corrections to the 
Einstein-Hilbert action of gravity at scales close to the 
4-dimcnsional Planck scale. However, the details of these 
corrections in a general time dependent background are 
less known. Thus, one would expect a generic action of 
type f{R,R^R^,R^ pa R^,ROR) 0,0,13. 

Of these examples, f(R) theories have received much 
attention due to their capability to mimic the late-time 
acceleration, see 0, 0, la 0, 0] including the solar sys- 
tem constraints On the other hand the nonlocal 
higher derivative corrections of the type R + RO l R 
yield an asymptotically free and ghost free nonpertur- 
bative action of gravity [lCj| . which has played an im- 
portant role in resolving the big bang singularity in the 
Friedmann-Robertson- Walker (FRW) universe. It also 
explains the observed temperature anisotropy in a non- 
inflationary bouncing universe setup [1 11 ] . Models of this 
type were also studied in [12], where it was shown that 
it is conformally equivalent to Einstein gravity coupled 
to two scalar fields. Also, models with nonlocal correc- 
tions but with negative powers of the d'Alambert opera- 
tor have been considered in [13| , where it was shown that 
such a theory may lead to the unification of early time 
inflation with late-time cosmic acceleration. 

However, in this paper we will not consider the anal- 
ysis for infinite, higher derivative nonlocal corrections of 
the type J2i RO l R, rather we will only concentrate on 
the i = l case, and the OR case. The complete nonper- 
turbative action will be dealt with separately in a future 
publication. 

The aim of the present paper is to study the low scale 
and long range behavior of a generic class of Lagrangian 
density, \f(R,R^ v R^ v ,OR), where we derive the scalar 
Newtonian potentials in a homogeneous and an isotropic 
expanding background such as in a FRW cosmology. The 
perturbations in the FRW background yield a Newto- 
nian potential for a matter distribution and therefore de- 
termine an effective Newtonian constant G e //. At long 
ranges the linear perturbation analysis differentiates Ein- 
stein's gravity with respect to any modification through 
the time evolution of the gravitational constant, see [lj], 



[lol ]. and this is one of the most important differences 
between Einstein and modified gravity theories . 

This difference can be tested by using the primordial 
nucleosynthesis (BBN) bounds on the gravitational con- 
stant, which are of the order of 10% [U S2, [H, Hi]. 
The BBN bounds are important due to the fact that the 
value of the gravitational constant determines the ex- 
pansion rate of the Universe and thus the relevant time 
scales for the production of light elements (H, He and 
Li), see [2^]. As a consequence, if we assume that the 
gravitational constant at the time of BBN is different 
from its value today, this means that the light element 
abundances will be different with respect to the standard 
BBN predictions. Even a weak time dependence, which 
gives no observable effects in Solar system experiments 
performed at the present epoch and at small scales, could 
give observable effects when translated over cosmological 
time scales. So, it will be interesting to analyze some 
special cases of our general Lagrangian and use the BBN 
bounds on the gravitational constant to place constraints 
on the parameters of these simple models. 

Furthermore, we will be applying the BBN constraints 
to study the OR corrections in the Einstein-Hilbert ac- 
tion. Previous studies of nonlocal action has concen- 
trated on formal aspects of the validity of effective field 
theory [25[ and particle creation [26[ . It should be noted 
that the effective Newton's constant in nonlocal gravity 
and its implications to cosmology, BBN and the Solar 
System have been also considered in |271| and more re- 
cently in models generalizing this in [281 ] . 

Here we consider the alteration of classical dynamics of 
the Universe due to the presence of the OR gravitational 
correction. 



1 The modifications in general relativity also affects structure for- 
mation |l6t UtI . [l8| and the predictions in the cosmic microwave 
background radiation through radiation-matter equality (l9, 20]. 
We will study various consequences to structure formation and 
Cosmic microwave background radiation in a separate publica- 
tion. 
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II. BACKGROUND EQUATIONS 

The action we will consider is 
1 



S 



d XyJ—g 



f(R,R^R^,OR)+L r , 



(2.1) 



where R is the Ricci scalar, R^ v the Ricci tensor, □ is the 
d'Alembert operator □ = g a @V a V ' p and L m the matter 
Lagrangian. We use the metric signature (— , +, +, +). 

Varying the action with respect to the metric g^ v we 
obtain the field equations as [29l| . [30l| : 



-(V(^) [V,)^)+^ 
df \ 



I? 



dUR 



[g^O - V„V„] (F + O^j + □ (f RR R^) + 

g^VcVp (f,RRR a0 ) - 2V a V (f iRR R% < } ) 



T, 



(2.2) 



where F = f£ and / RR 



Q R J,HH — d(R^R^) ■ 

defined the energy-momentum tensor as 



We have also 



T, 



2 8{y/=gL m ) 
~g 8g^ 



(2.3) 



and the parentheses next to indices mean symmetriza- 



tion, e.g. Aft) = \ (At 



A, 



Note that by using 



the field Eqs. (|2.2[) it is easy to see that it has the cor- 
rect limits, i.e. in the case when the Lagrangian is given 
by / = R + £OR then we get general relativity, as the 
OR term can be written as a total divergence. Also, for 
a conformally flat metric and a Lagrangian given bv[3~i| 
/ = .ff + £(3i? M „i? A " / — R 2 ) the field equations give general 
relativity at the background level, but not at the pertur- 
bations level as the metric is no longer conformally flat. 

In a flat FRW metric with a scale factor a(t), we obtain 
the zero-order (background) equations: 



/ ap/u2 , tt\ , nprj nff3 ( f „rj ni„_ d f \ , or tV tt2 ( d f 1 \ df 2 



2 - 3F ^ + H) + 3FH - 9F- ^2/ |Jul H - 2/, M - + Sfiff^ ^ + -f, RR ) - 6^H 

df , \,^vJi 9/ \ ^ 9/ £ ^ n .. T df 



-2^ - * - 3 (V;** + ^) ^ + ~ ^HrR - ^HfRR ~ 21H^ + H 

( ■ df \ df d 3 H d 4 9 { 

+ 6fRRH* - Af RR H - H Sf RR + 3^ - 24^/.,,. + - 4/, M _ + -g£ = p (2.5) 



where the dot ( ' ) denotes a derivative with respect to 
time, eg. f\ RR = and / iKH = g(R ^ Rllv) . 

III. PERTURBATION EQUATIONS 

We will consider the following perturbed metric with 
scalar metric perturbations $ and $ in a longitudinal 
gauge: 

ds 2 = -(1 + 2<$>)dt 2 + a{t) 2 {\ - 2^)8 ij dx i dx j (3.1) 

The energy-momentum tensor of the nonrelativistic mat- 
ter is decomposed as Tq = — (p m + Sp m ) and T° = 



—PmVm,a, where v m is a velocity potential. The Fourier 
transformed perturbation equations for the continuity 
equations are given by 

- PmV ™ k2 - 6 Pm - 3H6p m + 3*p,„ = (3.2) 

- a {Hp m v m + p rn u m ) = (3.3) 

Following the approach of Refs. [13], we use a sub- 
horizon approximation under which the leading terms 
correspond to those containing k 2 and Sp m . Terms that 
are of the form H 2 $ or $ are considered negligible rel- 
ative to terms like (fc 2 /a 2 )$ for modes well inside the 
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Hubble radius (fc 2 3> a 2 H 2 ). Under this approximation, 
the Fourier transformed perturbation equations, coming 
from the (/i, v) = (0, 0) and (1, 2) terms of the field Eqs. 
are given by 



Sp r , 



k 



— (SF-2F^) = 



(3.4) 



a 2 F 



dUR 



(* - 3*)/,jm = 



(3.5) 



While in general relativity in the case of a matter fluid 
with no anisotropic stress the two potentials $ and ^ are 
equal, as can be seen from Eq. (|3.5p . this is not the case 
for modified gravity theories as the gravity sector alone 
induces an anisotropic stress and creates the inequality 
of <3? and 4", see for example [HI, [33|. Next, we define the 
gauge invariant matter density perturbation 5 m as 



5 m = 



Spn 



3Hv 



(3.6) 



where 



v = av m (3.7) 
Under this approximation Eqs. (|3.2[) . (|3.3p and (|3.6p yield 

fc 2 $ 



2HS„ 



Next, we write 6F and S-^j^ as 







(3.8) 



(3.9) 



dUR 



dF 



SR 



+ 



dUR d{R llv R> iV ) 
d 2 f 



d4i- 

- R -j{R^in 



d{UR) 



AUR) 



(3.10) 



where 5R, under the subhorizon approximation, is given 
by 



5R~-2^(2V-<f>) 



while S(R^R^) ~ and 6 (OR) is given by 



(3.11) 



9 k 4 

5(DR) = (* - 2*) (3.12) 

a . 

Making these substitutions and using the subhorizon 
approximation in Eqs. (|3.4|) and (|3.5p we get 



k 2 fl 2 f k 8 



-2— - 2*) + f iRR ($ - *)) - 4 







dF k 6 
dUR^ 



($ - 2*) 
(3.13) 



i?($ - tf) + — ($ - 3*)/,«h - 2($ - 2*) 



fc 2 <9F fc 4 

~» + 2 



9 2 / k 6 



dDR a 4 d(DR) 



2 =0(3.14) 



The next step is to express $ and ^ in terms of 6 m . 
This can be done by solving the system of Eqs. (|3.13|) 
and (|3.14p for <f> and Doing so we find 



a 2 pS m F +M3/,flfl + 4F i 0-8^^+4 a( ^ )2 ^r 
" k 2 2(F - %f, RR ) F + £ (2/ fM + 3F fl ) - 6^£ + 3^£ 



(3.15) 



F+K (J tRR + 2Fj 



I dF fc 4 



o ,?7 , fc a 



fc 2 2(,F - ^f, RR ) F + f| (2/,«fl + 3*» - 6 



a_F fc 4 
ana a 4 



3- 



d 2 / fc" 



(3.16) 



9(D_R) 2 a 6 



From Eq. (|3.15p we can define a Poisson equation in constant G e tf. Doing so, we get the gravitational poten- 
the Fourier space and attribute the extra terms that ap- 
pear on the right-hand side to an effective gravitational 
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tial where G e // is defined as 



a 2 



$ = -47rG e// — S mPm (3.17) 
I 



ff ~RnF fe2 f ' P J- fc2 (Of -LIT? \ R 9F fc 4 in 3 2 / fc6 ^ J 



Note that the inclusion of the term RDR has a nega- 
tive contribution to G e / / . For certain choice of parame- 
ters it might be possible to make G e // vanishingly small, 
thereby modifying the Newtonian gravity on large tem- 
poral and spatial scales. 

Since the corrections from different forms of the mod- 
ifications, i.e. terms like R 2 , R^R^, OR etc enter with 
different powers of the k 2 it is interesting to check which 
hierarchies exist between the various coefficients in or- 
der for them to be equally important at some interesting 
scales. This can be very helpful to understand the rela- 
tive importance of the various modifications at different 
regimes. However, this is possible only for some simple 
cases and when the Lagrangian / is completely specified. 
In the general case it is not easy to tell whether a term 
of an arbitrary function, for example of RL~\R is more im- 
portant than some other term, as any of the derivatives 
of /, ie F, F : r etc may contain terms like DR. 

On the other hand by studying some simple cases, like 
the ones mentioned in the Examples section, we can draw 
some interesting conclusions. For example, as can be seen 
from Eq. 1|3.18[) for very small or very large scales ^ the 
terms containing DR are not as important as terms in- 
volving F and f,RR- However, on intermediate scales the 
DR terms can affect the behavior of G e // and actually 
enters with a negative sign which means that it may drive 
G e f f to zero or an unphysical singularity. 

Let us now study the Eq. (|3.8p of matter perturbations 



2H5 m - 4ttG, 



ef f Pm$m 



o 



(3.19) 



Note that the above expression will modify the large scale 
structure behavior on small scales as well as on large 
scales through higher order modifications. We will study 
these interesting possibilities in future publications. 



A. f{R) gravity 

As a first example we will consider f(R) theories, for 
which we have to set f(R, R^W, DR) = f(R). Then 
Eq. (|3~T5|) yields 



G 



eff 



1 l + 4^m 



(4.1) 



where 



RF 



m 



,R 



F 



being in agreement with the standard results from f(R) 
gravity [321 ]. 



f(R,R^R^) gravity 



A second example will be the Lagrangian 



f(R,R^R^,DR)=R 



UR^Rn n (4-2) 



In this case Eq. (|3.18[) gives 
G e //(a) 



8 ^ l-fE^^n^^)"- 1 

l + 2fEn M =0^(V)''- 1 



(4.3) 



If we keep only the first order term of the sum, corre- 
sponding to the Lagrangian R + ^R^R^, where £ is a 
constant, then G e // is 



IV. EXAMPLES 

In this section, we will consider several examples for 
the very general Lagrangian of the action (|2.ip in order 
to demonstrate how our results can be applied to a vast 
group of possible theories. 



G eff {a) 



1 + 3^ 



^ (1-^0(1 + 2^0 



(4.4) 



Next we will use the BBN constraints on the variation 
of the gravitational constant to constrain the parameter 
£. The effect of the variation of G e // can be constrained 
from BBN to be of the order of 10%, see, for example, 
Ref. [U, which gives ^fM^ = 1.09±°;f|. It is possible 
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FIG. 1: Plots of G e ff as a function of the scale factor a, given by Eq. (|4.4|l for example B (left) and by Eq. (|4.9[) for example 
C (right). The values of the parameter £ (in units of Mpc 2 and Mpc 6 respectively) used correspond to the best (blue line) and 
la values (green and red lines) allowed by the BBN bounds and are shown in each legend respectively. 



to use Eq. (|4.4[) to find analytically the best and the la 
values of £ according to BBN 



C. f(R, OR) gravity 



l BBN 



IBBN — 



Go 



± 



141 



Go 



In Fig. 1 (left) we show the plot of G e ff, given by Eq. 
(14. 3|) , for the values of the parameter £, which correspond 
to the central and la values allowed by the BBN bounds 
for a value of k = 0.002Mpc~ 1 . However, since the k- 
mode is actually unknown and can only be rather arbi- 
trarily chosen, we have also plotted the value of £ versus 
k in Fig. 2 and in Table 1 we show £ for various val- 
ues of the scale k. It is interesting to note that in this 
case there are actually two allowed values of £ by the 
BBN constraints, however only one is shown in Table 1 
for each k as the other results in completely unphysical 
behavior for G e //- Finally, we also consider the fc-mode 
corresponding to the horizon size at the BBN as the rel- 
evant scale. Since the horizon at the BBN is approxi- 
mately ~ 10~ 4 h~ 1 Mpc, see, for example, Ref. [3J], this 
corresponds to a scale ksBN ~ °f BN ~ 10~ 5 AIpc~ 1 and 
the corresponding constraints are snown in Table 1. 

It is possible to get more robust bounds on our models 
by considering local gravity constraints following the ap- 
proach of Ref. 32]. In this case we demand that strong 
modifications of gravity should not be observed on scales 
up to A/j ~ a/k, where in solar system experiments the 
scale Afe corresponds to a value around Xk — 1AU. There- 
fore, taking into consideration Eq. (|4.4p , we demand that 

fr|£| <C 1, which gives the following constraint 



As an example we will consider the case where the La- 

9( G q b o n ) 2 j grangian contains terms of the form DR. However, to 

keep the analysis simple we will consider only the first 

order term of such corrections and in this case the La- 
grangian will be given by 



(4.5) 



|£| «Af ~ IQ-^Mpc 



(4.6) 



While this is more robust than the ones found by us- 
ing the BBN constraint, the latter are not excluded as 
£bbn has a larger la error region, so the two constraints 
overlap with each other. 



f(R, R^R^, DR) = R + £DR 



(4.7) 



In this case the extra term DR can be rewritten as a 
total divergence and, as expected, does not contribute at 



TABLE I: The parameter £ using the BBN constrain for 
various values of the scale k. The first entry corresponds 
to the scale of the solar system experiments \k ~ 1AU or 
k so i ~ 2 10 11 Mpc' 1 , while the last (k ~ 10~ 5 Mpc' 1 ) corre- 
sponds to the scale of the horizon during the BBN. 



k (Mpc- 1 ) £ (Mpc 2 ) (case B) £ (Mpc 6 ) (case C) 



2 • 10 1 



If I « 10" 



If I « 10" 



1-10- 1 1.79*10- 17 ±^;*°:" 3.95 * 10- 48 ±^: 10 



10 -4S 



2-10- 3 4.48*10~ 14 ± 



14 I 1.04*10" 



9.46*10~ 14 



6.16 * 10 _ ^± 



38 I 2.15*10" 



1.00*10~ 37 



3-10" 4 1.99*10- 12 ± 4 : ^-i2 5.41 * io- 33 ±| : ?i;ig: 



1 ■ 10 1.79 * 10 ±3.7 8 „ 10 -9 3.95 * 10 ± 6 .4i*io- 
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FIG. 2: Plots of £ as a function of the scale k, given by Eq. (|4.5[1 for example B (left) and by Eq. (|4.10p for example C (right). 
The grey areas correspond to the la error bars. The scale that corresponds to the horizon at the BBN is Ubbn ~ 
and is situated outside of the range of the plots. 



all in the field equations. This can also be seen from the 
Friedmann Eq. (12. 4p , which in this special case simplifies 
to the usual Friedmann equation of Einstein gravity. 

The next most interesting case in this family of theories 
is the Lagrangian 



f(R, R^R^,OR) = R + £{OR) 2 



Then, G e // is given by 



G eff(a) 



1 1 



in l 



•6& 



(4.8) 



(4.9) 



All other cases involving terms (OR) n with n > 2 give 
complicated functions that also involve OR and thus are 
difficult to calculate. 

As in the previous example, it is possible to use Eq. 
(|4.9[) to find analytically the best and the la values of £ 
according to BBN 



l BBN 



,BBN 



1 I Gbbn i 

i+ G„ ) 



2fc 6 



(- 



Go 



(4.10) 



In Fig. 1 (right) we show the plot of G e ff, given by 
Eq. (|4.9j) . for the values of the parameter £, which corre- 
spond to the central and la values allowed by the BBN 
bounds for a value of k = 0.002Mpc~ 1 . However, since 
the fc-mode is actually unknown and can only be rather 
arbitrarily chosen, we have also plotted the value of £ 
versus k in Fig. 2 (right), and in Table 1 we show £ for 
various values of the scale k. Finally, as in the previ- 
ous case we will also consider the fc-mode corresponding 
to the horizon size at the BBN, and the corresponding 
constraints are shown in Table 1. 

Using the local gravity constraints for this example 
and taking into consideration eq. (|4.9|) . we demand that 
^■|f| <C 1, which gives the following constraint: 



|f | < XI ~ lQ-^Mpc 1 



(4.11) 



Again this is more robust than the ones found by us- 
ing the BBN constraint, the latter are not excluded as 
£bbn has a larger la error region, so the two constraints 
overlap with each other. 

Another very interesting case of this class of theories is 
to consider terms of the form ROR, instead of just OR. 
These terms correspond to the first order correction of a 
Lagrangian of the form R + J2^=o c n RO n R, which were 
shown in Ref. to give rise to a ghost and asymptoti- 
cally free theory of gravity. Thus, keeping only the first 
order correction the Lagrangian is 



f(R, Ry.uR^, OR) = R + £ROR 



and G e f / is given by 



G eff(a) 



l + fQR-8^ 



8tt(1 + (1 + £OR - 6^e) 



(4.12) 



(4.13) 



As can be seen by Eq. (|4. 13[) . G e ff also depends OR 
instead of just the scale factor a like in the previous cases. 
Unfortunately, we were unable to find cither an analytical 
solution, as the Friedmann Eq. (|2.4|) in this case is a very 
complex fourth order differential equation, or a numerical 
one as we do not have enough initial conditions. Thus, we 
were unable to provide a constraint for £ using the BBN 
bounds or plot G e ff as a function of the scale factor a. 



V. CONCLUSIONS 

Our analysis covers modified gravity models with a 
generic class of Lagrangian density with higher order 
and terms of the form ^f{R,R^ v R^ v ,OR). Using the 
fact that at long ranges the linear perturbation analy- 
sis differentiates Einsteins gravity with respect to any 
modification, through the time evolution of the gravita- 
tional constant, we derived the matter density perturba- 
tion equation and the effective gravitational "constant" 
G e f f for the action (|2.ip . 
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We also used the BBN bounds on the gravitational 
constant, which are of the order of 10%, in order to test 
the difference between Einstein and modified gravity the- 
ories. The reason why the BBN bounds can be used to 
test modified gravity theories is that the value of the 
gravitational constant determines the expansion rate of 
the Universe and thus the relevant time scales for the 
production of light elements (H, He and Li). This fact 
allowed us to test several cases of our general Lagrangian 
and constrain their parameters. Furthermore, we applied 
the BBN constraints to study the OR correction in the 
Einstein-Hubert action. 

However, the fact that the values we found for the 
parameter £ are actually larger than one would expect, 
it means that the energy scale at which these correction 
terms, e.g. R liV R! iV , are introduced is quite low. For 
instance one could write the corresponding term in the 
Lagrangian as j^R^R^" ', where ™ is the parameter £. 
Now, one would naively expect M to be of the order of 
Planck scale or even higher, but in our case the value of 
M is much smaller than that. 

This can be explained by the fact that presently the 
BBN bounds have quite a large error themselves, which 
means that the constraints we derived are not very 
strong. This can be seen by the fact that the error on 
the derived parameter £ is quite large, and this fact even 
allows for a zero value of £. Also, the primordial nu- 



cleosynthesis is quite a complex phenomenon and while 
its essence can be captured by a single data point, it 
is certain that a complete analysis, i.e. one that would 
also include the integration of the background equations 
from deep in the radiation era up to today and the use of 
the proper nuclear reaction rates, would most certainly 
provide stringent constraints. 

We have also implemented local gravity constraints, 
following the approach of Ref. *32j. As expected, the 
new constraints are more robust than the ones found by 
using the BBN constraint; however, the latter are not 
excluded as {,bbn has a larger la error region, so the 
two constraints overlap with each other. 
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